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ABSTRACT 

The M dwarf V374 Peg (M4) is believed to lie near the theoretical mass threshold 
for fully convective interiors. Its rapid rotation (P = 0.44 days) along with its intense 
magnetic field point toward magneto-centrifugal acceleration of a coronal wind. In this 
work, we investigate the structure of the coronal wind of V374 Peg by means of three- 
dimensional magnetohydrodynamical (MHD) numerical simulations. For the first time, 
an observationally derived surface magnetic field map is implemented in MHD models 
of stellar winds for a low mass star. By self-consistently taking into consideration the 
interaction of the outflowing wind with the magnetic field and vice versa, we show 
that the wind of V374 Peg deviates greatly from a low-velocity, low-mass-loss rate 
solar-type wind. We have found general scaling relations for the terminal velocities, 
mass- loss rates, and spin-down times of highly magnetized M dwarfs. In particular, for 
V374 Peg, our models show that terminal velocities across a range of stellar latitudes 

reach Uqq ~ (1500 - 2300)n 12 1/2 km s~\ where ni2 is the coronal wind base density in 

units of 10 12 cm -3 , while the mass-loss rates are about M ^ 4 x lO -10 ?^ 2 M yr _1 . 
We also evaluate the angular- momentum loss of V374 Peg, which presents a rotational 

— 1/2 

braking timescale r ~ 28n 12 Myr. Compared to observationally derived values from 
period distributions of stars in open clusters, this suggests that V374 Peg may have low 
coronal base densities (< 10 11 cm -3 ). We show that the wind ram pressure of V374 Peg 
is about 5 orders of magnitude larger than for the solar wind. Never the less, a small 
planetary magnetic field intensity (~ 0.1 G) is able to shield a planet orbiting at 1 AU 
against the erosive effects of the stellar wind. However, planets orbiting inside the 
habitable zone of V374 Peg, where the wind ram pressure is higher, might be facing 
a more significant atmospheric erosion. In that case, higher planetary magnetic fields 
of, at least, about half the magnetic field intensity of Jupiter, are required to protect 
the planet's atmosphere. 

Key words: MHD - methods: numerical - stars: individual (V374 Peg) - stars: 
low-mass - stars: magnetic field - stars: winds, outflows 



1 INTRODUCTION 

M dwarf (dM) stars comprise the most abundant group of 
stars in the Galaxy. Due to their small masses (~ 0.06 - 
0.8 M , Sca lo et aLlboO/t ). they have long main-sequence 
lifetimes, which exceed the Hubble time. They are low- 
luminosity objects and because of that, their habitable zone, 
where liqui d water may be foun d, is believed to lie very close 
to the star ([Kasting et alJll9 93). This makes dM stars prime 
targets in searches for terrestrial habitable planets. 

From the point of view of stellar evolutionary theory, 
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dM stars lie on the boundary between fully convective stars 
(spectral types later than ^M4) and stars with a radiative 
core (spectral types earlier than ^M4). The latter are be- 
lieved to host magnetic fields generated by solar-type dy- 
namos, where the layer between the radiative core and the 
convective outer shell (the tachocline) produces rotational 
shear, which, in conjunction with convective motions, is able 
to amplify and maintain surface magnetic fields. The genera- 
tion of magnetic fields in fully convective stars is still a mat- 
ter of debate though. Because such stars lack a tachocline, 
the dynamo mechanism should be different from a solar- 
type dynamo. Also, their nearly solid-body rotation indi- 
cates that the magnetic fields should cause a quenching of 
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the surface differential rotation, an ingredient believed to be 
essential in field generation. 

Although details of the dynamo mechanism may still 
be unknown, the change in the stellar internal condi- 
tions appears to be reflect in the shift in the magnetic 
field large-scale topology around spectral type M4. The 
observation of surface magnetic field distributions sug- 
gests that early-dM stars host weak large-scale fields with 
dominantly toroidal and non- axisymmetric poloidal con- 
figurations ([Donati et alj [2008), while mid-dM stars host 
strong, mainly ax i-symmetric large-scale poloidal fields 
([Morin et al.ll2008al ). 

This topology shift should produce a corresponding 
change in the coronal structure of the star, and, in partic- 
ular, in coronal tracers such as X-ray and radio emissions. 
However, both the X-ray luminosity Lx as well as the quies- 
cent radio luminosity L ra d do not present sharp transitions 
across the M4 boundary. L ra d, for instance, remains rela- 
tively the same for all the spectral subtypes of dM stars, 
while Lx is roughly proportional to the bolometric lumi- 
nosity Lboi, presenting a s aturation limit of Lx jLhoi ~ 10~ 3 
out t o spectral types M7 (Delfo sse et al.l Il998l ; IBerger et al.l 
2010). As a consequence, a correlation between radio and 
X-ray emissions is observed for almost all the spectral sub- 
types of dM stars, with a s harp transition b e ing seen only 
around spectral ty pe M7 ( Giidel et all Il993l : iBergerl 12006 : 
IBerger et al.ll201Q[ ). Because of such a correlation, it is often 
suggested that the origins of both X-ray emission, which is 
believed to trace the hot plasma contained in magnetic loop 
structures, and radio emission, which is believed to be pro- 
duced by gyro- synchrotron radiation, are common and most 
probably due to the coronal magnetic fields. 

Regarding rotation, the evolution of dM stars can be 
inferred from obse r vations of open clusters at different ages 
(llrwin et al.1 120061: [Scholz EislofTell 120071: iHartman et al.l 
2009 : ICollier Cameron et al.ll2009l : iMeibom et al.1 120091 ). In 
young (< 700 Myr) open clusters, dM stars still present 
high rotation rates, which suggests that angular momentum 
losses at the early main-se quence phase are negligible for 
them (|lrwin &; Bouvierll2009h . However, as the cluster ages 
(> 700 Myr), the number of rapidly rotating dM stars de- 
creases, implying that there should exist a mechanism of an- 
gular momentum remo val that acts on time-s cales of a few 
hundred million years (|Scholz fe Eisloffel 20071 ). For solar- 
like main sequence stars, the magnetised stellar wind is be- 
lieved to spin down the star by carrying away stellar angular 
momentum. It has been observationally established that the 
angular velocity rate Qo for solar- like stars va ries as a func- 
tion of age t as Oo oc t -1 / 2 (ISkumanichlll972l ). However, it 
seems that the empirical ISkumanichl 's law is not valid for 
low-mass stars, suggesting that a solar- type wind (i.e., with 
low velocities and mass-loss rates) cannot reproduce the ro- 
tational evolution of fully- convective stars. 

The existence of hot coronae, rapid rotation, and high 
levels of magnetic activity in dM stars suggests the presence 
of winds with an enhanced mass loss as compared to the 
solar wind. However, the low-density, optically thin winds of 
these stars prevents the observation of traditional mass-loss 
signatures, such as P Cygni profiles. The still unobserved 
high mass-loss rates from dM stars could be able to disperse 
debris discs, explaining why discs around dM stars older 
than > 10 Myr are scarcely found ([Plavchan et al.l 12005). 



Estimates of mass-loss rates from dM stars vary 
considerably. It has been suggested that the coronal 
winds of dM stars, despite of being very tenuous, pos- 
sess mass-loss rates (M) that can considerably exceed 
the solar value ( M ~ 2 x 10" 14 M yr" 1 ) by fac- 
tors of 10 to 10 4 dMullan et al.lll992l: iBa dalv an &; Livshits 



Il992l: iLimfe White! Il996l : Ivan den Oord Dovld Il997l : 
IWargelin fc Drake 200ll ). Based on radio observations from 
IRAS, VLA and JCMT data. lMullan et al.1 (|l992h estimated 
mass-loss rates for active dM stars (in particular, YZ CMi, 
Gl 644, Gl 873) of a fe w times IP' 10 M yr" 1 . This re- 
sult was contradicted bv lLim fc White] (l996), who claimed 
that the active M dwarfs YZ CMi (dM4.5e) and AD Leo 
(d M3.5e) should present M < 10 ~ 13 M yr -1 , and later on 
by Ivan den Oord h Dovld (|l997l ). who suggested that the 
mass-loss rates of dM stars should be at maximum acou- 
ple of M < 10" 12 M yr" 1 . Ivan den Oord fe Dovld (|l997l ) 
pointed out that both YZ CMi and AD Leo are consid- 
ered young stars (ages of about 0.5 - 1 Gyr), and that 
at their early life, M should have been about one order 
of magnitude larger: M ~ (1.1 - 3.6) x 10" 11 M yr _1 . 
Even for the closest star to the Solar System, Proxima Cen- 
tauri (dM 5.5e), estimates of mass- lo ss rates are rather con- 
troversial. Warge lin fe Drakd ([2002) suggest an up per limit 
of M < 3 x 10" 13 M yr _1 and iLim et all (jl996l ) suggest 
M < 7 x 10" 12 (^oo/300 km s" 1 ) M yr" 1 (where is 
the wind terminal velocity), although IWood et al.1 (2001) 
claim an upper limit of M < 4 x 10 -15 M yr _1 , 5 times 
below the value of the solar wind mas s-loss rate. Accord- 
ing tolz endeias, Segura, fe Ragal (|2010h . if Proxima Cen is 
about 1-4 Gyr old, the "young" Proxima Cen should 
have had considerably larger mass-loss rates M ~ (4 - 
64) x 10" 11 M yr _1 . 



In this work, we investigate the coronal wind of a 
specific fully-convective dM star, V374 Peg, for which 
observed sur f ace magnetic map s have been acquired 
(jDonati et al.1 120061 : iMorin et all l2008bl ). For this, we 
use three-dimensional (3D) magnetohydrodynamics (MHD) 
simulations based on our previous m odels developed for 
solar-like stars (IVidotto et al.1 120 09a) and weak- lined T 
Tauri stars (| Vidotto et al.ll2009bl . l201oh . For the first time, 
an observationally derived surface magnetic field map is im- 
plemented in MHD models of stellar winds for a low-mass 
star. V374 Peg is a suitable case for modelling as a first step, 
because its surface magnetic distribution is close to poten- 
tial, which implies that the adopted boundary conditions 
match the observed map closely. We self-consistently take 
into consideration the interaction of the outflowing wind 
with the magnetic field and vice-versa. Hence, from the in- 
terplay between magnetic forces and wind forces, we are able 
to determine the configuration of the magnetic field and the 
structure of the coronal winds. 



This paper is organised as follows. In Section [2] we 
present the details of the stellar wind model we use to de- 
scribe the wind of V374 Peg. Section [3]presents the results of 
our simulations. In Section 21 we present the discussions of 
our work, including implications of the wind on habit ability 
of planets. Section [5] presents the conclusion of our work. 
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2 THE NUMERICAL STELLAR WIND MODEL 

Observations of the solar corona during eclipses have re- 
vealed the presence of closed magnetic loops, extending 
out to several solar radii. Beside the closed loops, re- 
gions of open magnetic field lines form coronal holes, from 
where the fast solar wind originates. Based on the observed 
solar configuration, bi-component coronal models, where 
closed and open magnetic field line regions coexist, have 
been used t o explain the coronal st r ucture of the Sun and 
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In these models, closed field structures form a belt around 
the equator of the star, while open field structures give rise 
to the fast solar /stellar wind. We adopt such wind models 
to simulate the yet unobserved wind of V374 Peg. 

To perform the simulations, we use the 3D MHD nu- 
meri cal code BATS-R-U S developed at University of Michi- 
gan flPowell et al.lll999h . BATS-R-US solves the MHD equa- 
tions in the finite- volume scheme. It has currently several 
approximate Riemann solvers available, although originally 
used the upwind Roe- type scheme. The temporal integration 
is performed using either an explicit, implicit or combined 
explicit-implicit time stepping scheme (details are found in 
iToth et alJl2006h . We perform our calculation in an explicit 
scheme, where a two-level Runge-Kutta scheme is used. 
Here, we summarise th e characteristics of o ur model. More 
details can be found in lVidotto et al.l (|2010l ). 

BATS-R-US solves the ideal MHD equations, that in 
the conservative form are given by 
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where p is the mass density, u the plasma velocity, B the 
magnetic field, p the gas pressure, g the gravitational accel- 
eration due to the central body, and e is the total energy 
density given by 



pu 



1 + 8tt" 



(5) 



We consider an ideal gas, so p = nkBT, where Ub is the 
Boltzmann constant, T is the temperature, n — p/(pm p ) is 
the wind density, pm p is the mean mass of the particle, and 
7 is the heating parameter (polytropic index). In practice, 
7 represents the energy balance of the wind, for which the 
details are unknown. In our simulations, we adopt 7 = 1.2 
and p — 0.5. 

BATS-R-US uses a Cartesian computational domain 
that is block-based, where blocks can be either created or 
destroyed as the computation proceeds. To do that, an adap- 
tive mesh refinement (AMR) technique is used in the MHD 
algorithm of BATS-R-US. Although we opted not to mod- 
ify our grid resolution during the numerical computation, 
we used the AMR capabilities of BATS-R-US to construct 



a grid for the simulations of the wind of V374 Peg that is 
more refined around the star, which is placed at the centre 
of the Cartesian grid. The grid axes x, y, and z extend from 
—20 i?* to 20 i?*, where R* is the stellar radius. After apply- 
ing 8 levels of refinement, the more refined zone of the grid 
(for r < 5 i?*) has cells with a resolution of 7.8 x 10 -2 i?*. 

The inner boundary of the system is the base of the 
wind at radius r = R*, where fixed boundary conditions 
are adopted. The outer boundaries at the edges of the grid 
have outflow conditions, i.e., a zero gradient is set to all the 
primary variables (u, B, p, and p). 

V374 Peg has mass M* = 0.28 M©, radius R* = 
0.34 Rq and is rotating with negligible differential rotation 
i.e., as a solid body ) with a period of rotation Po = 0.44 d 
Morin et al.ll2008bh . We consider that its axis of rotation 
lies in the ^-direction. The simulations are initialised with a 
ID hydrodynamical wind for a fully ionised plasma of hydro- 
gen. Immersed in this wind we consider a magnetic field an- 
chored on the stellar surface that has either a dipolar geom- 
etry ( N2.ip or a geometry derived from extrapolations from 
observed surface magnetic maps f ^2.2[) using the potential- 
field source surface (PFSS) method. The MHD solution is 
evolved in time from the initial magnetic field configuration 
to a fully self- consistent solution. We do not adopt fixed 
topologies for either the magnetic field or for the wind, as 
both the wind and magnetic field lines are allowed to inter- 
act with each other. 



2.1 Dipolar Surface Field 

Because V374 Peg presents a surface magnetic field distri- 
bution that is mainly dipolar, a dipolar field geometry was 
initially used for our simula tions. This allows com p arisons 
with the previous works of IVidotto et ail (|2009al lbl l2010h 
and with a realistic magnetic field distribution at the stellar 
surface. In such cases, the initial magnetic field configuration 
is described in spherical coordinates (r, ip) by 



BqR* 



cos 6, - sin#, 



(6) 



To be compatible to the maximum intensity derived in the 
surface maps, we adopted a maximum intensity of the mag- 
netic field of Bo = 1660 G, which is evaluated at the 
magnetic pole. Here, we assume that the dipole moment is 
aligned with the rotation axis of the star. Figure [TJi shows 
the initial configuration of the magnetic field lines. 



2.2 Observed Surface Field 

In the remaining cases, we implemented a more realistic dis- 
tribution of magnetic fields at the surface of the star. This is 
the first time an observed magnetogram has been included 
in a MHD wind model for a low- mass star. Magnet i c map s 
for V374 Peg have been presented in iDonati et al.1 (|2QQ6h : 
I Morin et al.l ( 2008bh , and are derived from data acquired in 
2005 Aug, 2005 Sep, and 2006 Aug. As a first step, we have 
opted to implement the maps derived from the particular 
set of observations made in 2005 Aug. Figure [2] shows the 
radial, meridional and azimuthal components of the surface 
magnetic field of V374 Peg obtained through extrapolations 
of the surface map using the PFSS technique, as described 
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(a) Dipole: Initial configuration 



(b) Surface map: Initial configuration 




(c) Dipole: Final configuration 





x(R.) 

(d) Surface map: Final configuration 






x(R.) 

Figure 1. (a) Initial configuration of magnetic field lines for a dipole and (b) for the extrapolation of the surface map using the potential- 
field source surface (PFSS) technique, where the dashed circumference represents the position of the source surface, (c) Final configuration 
of the of magnetic field lines after the self-consistent interaction with the stellar wind considering a dipolar surface distribution of magnetic 
fields and (d) considering the observationally derived magnetic map. Note that panels (b) and (d) are shown for different stellar rotation 
phases. 



in ljardine et al.1 (|2002h . In the PFSS model, the stellar wind 
plasma is not included directly, but its effects on the mag- 
netic field (and vice- versa) are incorporated through the in- 
clusion of the source surface. Such a surface, for instance, al- 
ters the number of open magnetic field lines, through where 
a stellar wind could escape. PFSS methods are usually crit- 
icised because their basic assumptions (the magnetic field is 
a potential field and the source surface is spherical) may not 



always be met. However, the advantage of the PFSS method 
over the MHD models relies on its simplicity: it is simpler to 
implement and requires much less computer resources than 
MHD models. In our model, we use the magnetic field con- 
figuration derived by the PFSS method as initial condition 
and boundary condition at the surface of the star. We note 
that the surface maps presented in Figure [2] matches the 
observed maps closely. 



Powerful Winds from Low-Mass Stars: V374 Peg 5 



Radial 



135 



300 360 




180 



120 180 240 

Longitude (degrees) 

Azimuthal 



300 360 



2 45 



a 90 



T 135 




120 180 240 

Longitude (degrees) 



360 



Figure 2. The three components of the surface magnetic field 
of V374 Peg obtained through extrapolations of the surfac e map 
using t he PFSS technique. Maps were derived by iDonati et al.l 
( 2006): [Morin et"aTI (l2008bh an d extrapolations were done as de- 
scribed inTja3In^^niiQ02) • 



A meridional cut of the initial configuration of the mag- 
netic field lines is presented in Figure [TJ). The source sur- 
face (dashed circumference) is chosen to lie at rss = 5 R*, 
where beyond that, the magnetic field is considered to be 
purely radial. We can see that the extrapolated field is in- 
deed mainly dipolar, slightly inclined with respect to the 
rotational axis of the star (z-axis). We note that the surface 
of the star that occupies co- latitudes > 120° is never in view 
as the star rotates and so the magnetic field there can not 
be reconstructed reliably. 



3 RESULTS 

Our simulations require a set of input parameters for the 
wind. Unfortunately, some of them are poorly constrained by 
observations. For V374 Peg, the magnetic field is the better- 
constrained paramete r. We have, therefor e, implemented in 
our previous models (|Vidotto et al.| [2010) sur face magnetic 
maps derived by data acquired in 2005 Aug (|Donati et al.l 
2006). These observations show that V374 Peg hosts an in- 
tense, mainly axi-symmetrical dipolar magnetic field, with 



maximum intensity of about 1660 G, i.e., 3 orders of mag- 
nitude larger than the Sun. 

The wind temperature and density are less constrained 
for V374 Peg. We, therefore, adopt values representative 
of dM stars. dM stars are believed to host coronae with 
a high-temperature plasma ~ 10 7 K in conjunction with 
a low-temperature on e 2 - 3 x 10 6 K ([Schmitt et al.| [l990; 
iGiampapa et al. 1996). In our models, we adopt a temper- 
ature at the base of the coronal wind of To = 2 x 10 6 K 
or 10 7 K. These coronal temperatures are about the same 
order of magnitude as the solar coronal temperatures of 
1.56 x 10 6 K. 

Coronal densities inferred from X-ray observations of 
dM stars sug gest densities rangin g from 10 10 cm -3 to 5 x 
10 12 cm" 3 (jNess et al.ll200ll2004l ). Therefore, we adopt, at 
the base of the coronal wind, densities in the range 10 10 - 
10 12 cm -3 . Compared to the solar coronal density of about 
~ 2 x 10 8 cm -3 , coronal densities inferred for dM stars are 
about 2-4 orders of magnitude larger than for the solar 
corona. 

The density, along with the magnetic field, are key pa- 
rameters in defining the magnetic f ield configuration of the 
stellar wind and its velocity profile ([Vidotto et al.ll2009al lbh. 
Together, they define the plasma- /3, defined by the ratio of 
thermal to magnetic energy densities. Therefore, at the base 
of the coronal wind of V374 Peg, 



_ no kB Tp 



■ 2.5 x 10 5 n 10 T 6 , 



(7) 



where the index "0" means the variable is evaluated at 
the base of the coronal wind, mo = no/(10 10 cm -3 ) and 
T Q = T /(2 x 10 6 K). For mo = 1 and T 6 = 1, /3 is 
about 5 orders of m agnitude smaller than for the solar wind 
([Pneuman &; Kopplll97ll ) . This implies that the winds of dM 
stars are highly magnetised and, therefore, are expected to 
differ from solar- type winds. 

Table [1] presents the parameters adopted for the set of 
simulations we performed. The suffixes 'Dip' and 'Map' at 
the case name (first column) stand for the distribution of 
the surface magnetic field adopted: dipolar or from a mag- 
netic map, respectively. The dipolar cases form the Set 1 of 
simulations while the cases with magnetograms form the Set 
2. 

We were able to find a MHD solution for the wind for 
all the simulations we ran, showing that it is possible to de- 
velop coronal wind models with a realistic distribution of 
magnetic field. In general, MHD wind models are studied 
under the assumption of simplistic magnetic field configura- 
tions, especially when in pursuit of an analytical solution. 
Therefore, the study of a magnetised coronal wind where an 
observed magnetic field distribution is considered has long 
been awaited. Our work also sheds some light on the yet 
unobserved winds from dM stars. 

Our results show that cases where all parameters are 
the same, except for the surface magnetic field configura- 
tion (e.g., IDip/lMap, or 2Dip/2Map, and so on) present 
similar general wind characteristics and magnetic field con- 
figurations. The wind velocities are approximately the same 
for cases 3Dip and 3Map, although the details of the velocity 
profiles may differ. Overall, our solutions differ considerably 
from the solar wind solution, where a low- velocity wind (ter- 
minal velocities of Uqo,q ~ 400 - 800 km s" 1 ) with low mass- 
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Table 1. Adopted parameters for the sets of simulations. The columns are, respectively: the case name, the density no and temperature 
To at the base of the coronal wind (r = i?*), the plasma- (3 parameter evaluated at R* [Eq. [7], the mass loss rate M, the angular 
momentum loss rate J, and the time-scale for rotational braking r. The case names are labelled either "Dip" or "Map" to indicate that 
the distribution of the surface magnetic field is, respectively, dipolar or from an observed magnetogram. 



Case 


no 


To 


/3o 


M 


j 


T 




(cm- 3 ) 


(MK) 


(io- n 


M yr- 1 ) 


(10 33 erg s- 1 ) 


(Myr) 


Set 1: Dipole 


IDip 


10 10 


2 


2.52 x 10- 5 


3.4 


2.3 


280 


2Dip 


10 11 


2 


2.52 x 10- 4 


12 


7.7 


84 


3Dip 


10 12 


2 


2.52 x 10-3 


41 


26 


24 


4Dip 


10 11 


10 


1.26 x 10- 3 


25 


8.3 


78 


Set 2: Map 


IMap 


10 10 


2 


2.52 x 10~ 5 


4.2 


3.4 


180 


2Map 


10 11 


2 


2.52 x 10~ 4 


14 


7.6 


84 


3Map 


10 12 


2 


2.52 x lO -3 


50 


32 


17 


4Map 


10 11 


10 


1.26 x lO" 3 


26 


9.1 


48 



loss rate (M© ~2x 10~ 14 M yr -1 ) is found. Terminal ve- 
locities of our models are around Uoo ~ 1500 - 2300 km s _1 
for cases where no = 10 12 cm -3 , increasing as no decreases, 
while the mass- loss rates are about M « 4 x 10 -10 M© yr -1 
increasing as no increases. A comparison of terminal ve- 
locities derived from simpler wind models is presented in 
Appendix [Aj We note that, based on more simplistic wind 
models, such as IWeber fe Davisl <|l967h , in the fast magnetic 
rotator limit, wind terminal velocities of ~ 3320 km s _1 are 
expected for a wind mass- loss rate of about 10 -11 M© yr -1 . 
Figures [Tt and[TJi present the final configuration of the mag- 
netic field lines, where we note the similarities of the field 
topology for cases with a dipolar distribution of surface mag- 
netic flux (Fig. [I]:) and with a distribution given by observed 
magnetic maps (Fig. [TJi). As a result of the rotation of the 
star, the apparent positions of the magnetic poles in the 
initial (Fig. [T]d) and in the final magnetic field configura- 
tions (Fig. [TJi) are different because the plots represent two 
differ ent stellar phases of rotation. 

M orm et al.) (|2010h suggested that the magnetic field 
structure of very low-mass stars may switch from a strong, 
simple field configuration to a weak, complex field over time. 
If this is the case for V374 Peg and its magnetic field is in 
such a transient state, its wind characteristics as well as 
mass-loss and angular momentum-loss rates are expected to 
change along with its magnetic field structure. 

In the next sections, we present more details of our re- 
sults divided in two sets defined by the configuration of the 
surface magnetic distribution adopted. 



3.1 Set 1: Dipolar Field 

Set 1 considers a magnetic surface distribution of a dipolar 
field. By comparing cases IDip, 2Dip, and 3Dip, i.e., where 
only the base coronal density no was varied, we found that 
the poloidal velocity of the wind scales approximately as 



ui oc — for a given Bo. 
n 



(8) 



This qualitatively agrees with previous results 
(Vid otto et all [2009a), where it was found that an in- 



crease in the density leads to winds with lower velocities. 
Simulations IDip, 2Dip, and 3Dip are in a very low-/? 
regime, where magnetic effects completely override thermal 
and kinematic effects of the wind. The equation of motion 
in steady state reduces to 



p(u • V)u = -Vp + pg + 



j x B 



(9) 



where j = c(V xB)/ (4ty) is the current density. For simplic- 
ity, let us take the radial component of previous equation 
and integrate it. Neglecting terms with p and pg relative 
to the much larger magnetic term, the previous equation is 
rewritten as 



/ 



\j x B|. 

cp 



-dr. 



(10) 



This shows that in a domain where p <C B 2 (very low- 
f3 regime) and gravitational forces can be neglected, the 
squared velocity of the wind is inversely proportional to the 
density [Eq. ©]. 

At the Alfven surface, the poloidal wind velocity equals 
the local Alfven velocity 



2 B p a 

U p ,A - 



(11) 



As u 2 A oc 1/pA for cases IDip, 2Dip, and 3Dip [Eq. JS}], 
the location of the Alfven surface will only depend on the 
intensity B p ,a of the magnetic field at the Alfven surface. As 
these cases possess the same magnetic field at the base of the 
coronal wind, and because the configuration of the magnetic 
field lines is very similar, the magnetic field at the Alfven 
surface should be approximately the same. This implies that, 
despite the increase in the base density, the location of the 
Alfven surface is similar for cases IDip, 2Dip, and 3Dip . 
This is different to what was found in I Vidotto et aT1(|2009ah . 
where an increase in the coronal density modified the loca- 
tion of the Alfven surface. This Alfven surface saturation 
happens because simulations IDip, 2Dip, and 3Dip are in a 
very \ow-/3 regime. Therefore, Eq. (|SJ) should be treated with 
caution, as under different ft regimes (for example, when it 
approaches /3q ~ 1), it becomes invalid. Figure [3^ presents 
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Figure 3. Meridional cut of scaled poloidal wind velocity u p = u p / \/n\i profile (a) for cases with a dipolar field surface distribution 
IDip, 2Dip, and 3Dip, and (b) 4Dip, (c) for cases where a magnetic surface map was used IMap, 2Map, and 3Map, and (d) 4Map. Black 
lines represent the magnetic field configuration, and white line is the location of the Alfven surface. 



the scaled wind velocity profile u p for cases IDip, 2Dip, and 
3Dip, 



1/2 



(12) 



We note that the wind ter- 

-1/2 „-l 



where nn — no/(10 12 cm 

minal velocity is Uoq ~ (1500 - 2300) n^ 1 1 z km s _1 , where 
the range of velocities refers to different wind latitudes (low- 
wind velocity near the equator, high- wind velocity around 
the poles). 

Equation also implies that the mass-loss rate of the 
wind (M oc pu r ) should scale as 



M oc pu r oc n 



1/2 



(13) 



which means that, despite the fact that the wind velocity 
of case 3Dip is 10 times smaller than case IDip [Eq. 0], 
its mass-loss rate is one order of magnitude larger than for 
case IDip [Eq. (| 13)) ] . This has implications on the efficiency 
of angular momentum loss, as will be shown in fc|4.11 The 



mass-loss rates for cases IDip, 2Dip, and 3Dip are M w 
4 x 10- 10 n^ 2 M yr- 1 . 

Figure [3J) shows the scaled poloidal velocity profile u p 
for case 4Dip. This case considers a different temperature at 
the base of the corona (10 7 K as opposed to 2 x 10 6 K), and, 
because of that, has a larger f3o (Table [T]). For this case, we 
did not find an analytical expression relating velocity and 
temperature. The Alfven surface location and configuration 
of magnetic field lines are similar to the other dipolar cases, 
but the scaled wind velocity u p is smaller than the previous 
cases: i^oo ~ (800 - lSOO)™^ 2 km s _1 . The lower velocity 
observed in case 4Dip happens because of its higher /3o. The 
mass- loss rate for case 4Dip is M ~ 2.5 x 10 -10 M yr _1 . 



3.2 Set 2: Observationally Derived Magnetic Map 

The results of the simulations where the magnetic field dis- 
tribution at the surface of the star is set according to ob- 
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served magnetic maps is presented in a similar fashion as 
done in 33.11 the profile of the scaled velocity u p [Eq. (p"2|) ] 
for cases IMap, 2Map, and 3Map is shown in Figure [3k, 
while for case 4Map it is presented in Figure [3ji. Because 
the magnetic field in the lower hemisphere of the star is not 
reliably reconstructed (co- latitudes > 120° of the surface of 
the star are not observed), a high- velocity wind develops 
there. Although this feature is local and does not affect the 
remaining parts of the grid other than radially away from 
the stellar surface, it is an artifact of our method and should 
not be taken into consideration (e.g., see the accumulation 
of magnetic field lines near the low- hemisphere of the star 
in Figure [TJ)) . The wind velocities of Set 2 are consistent 
with the wind velocities of Set 1, shown in 33.11 This can 
be seen by comparing Fig. [3^ and [3t, and [3J) and [3ji. We 
found similar values of wind terminal velocities (i^oo ~ [1300 
- 2100]n7 2 1/2 km s -1 for cases IMap, 2Map, and 3Map, 
and i^oo ~ [850 - lGOOjn^ 2 km s _1 for case 4Map) and 
mass- loss rates (M « 4 x lO -10 ^ 2 M vr" 1 for cases 
IMap, 2Map, and 3Map, and M - 2.6 x 10" 10 M yr" 1 
for case 4Map). This similarity occurs because the observed 
magnetic map does not deviate significantly from the dipo- 
lar field adopted in 33.11 Future works will consider other 
stars with more complex surface magnetic field distributions, 
where more significant deviations are expected. 

Tests with different positions for the source surface 
showed that final configuration of the wind does not differ 
from what is presented in this paper. 



4 DISCUSSION 

4.1 Constraining Model Parameters: Angular 
Momentum Loss 

One way to restrict the parameters of our model is by 
comparing the derived values for Uoo or M with ob- 
served values. However, traditional measurements of mass- 
loss rates or wind terminal velocities have not yet been 
precisely obtained for dM stars. Estimates of mass-loss 
rates from dM stars, for instance, are rather controver- 



loss rate (e.g., iMullan et a 


[J 19921: Badalvan & Livshitd 


1992: iLim fc White! Il996l: 


van den Oord & Dovld Il997l: 



IWargelin &; Drake! 1200 lh, but pred ictions for sub- solar val- 
ues exist as well ([Wood et al.ll200lh . Mass- loss rates derived 
from our wind models are M « 4 x 10" 10 n^ 2 M vr -1 
for cases IDip/lMap, 2Dip/2Map, and 3Dip/3Map, and 
M - 2.5 x 10" 10 M vr -1 for case 4Dip/4Map. Such val- 
ues are in accordance with several estimates. Unfortunately, 
they do not allow us to constrain the input parameters of 
our model, in particular, the density at the coronal base. 

Another possible way to restrict our model parame- 
ters is to compare the time-scale for rotational braking r 
with the age t of V374 Peg. Because V374 Peg is a fast 
rotator, we expect that r ^> t. However, there exist only 
loose age estimate s for V 374 Peg. iDelfosse et aTl fl998) and 
iMohantv &; Basri (|2003h estimated the kinematic age of 
V374 Peg as "young disc", which roughly corresponds to 
t < 3 G yr. Using the calibration of Lx versus age pro- 
vided by ICatalan etHI (l2009h. for the observed value of 
L x = 10 28 - 44 erg s _1 (jPelfosse et al.lll998h . it results in 



an age estimate o f t ~ 0.7 - 0.8 Gyr. The calibration from 
West et alJ (|2009h for M4 dwarfs results in the considerably 
older age of t ~ 5.4 Gyr, using Ln a and Lboi as given by 
IDelfosse et al.l (|l998l ). Because age estimates are still loose 
and imprecise, its comparison with r is not able to restrict 
our model parameters. 

Observations of the rotation evolution of dM stars in 
open clusters at different ages, on the other hand, provide a 
way to constrain the time-scale r for the angular-momentum 
loss. It has been suggest ed that r ^ 200 M yr or, mostly likely 
400 - 800 Myr, (jScholz fc Eisloffelll2007l ) for dM stars. 

Angular momentum of the star is carried away by the 
stellar wind. Because in some of our simulations there is no 
axi-symmetry, the torque J on the star has x, y and z com- 
ponents. Here, we are interested only on the z-component, 
as it is the only one responsible for the rotational brak- 
ing (because the angular velocity of the star points in the 
z-direction). The ^-co mponent of the angul ar momentum 
carried by the wind is (jMestel fe Sellevlll970h 



J = 



+ 



/ 

Jv A 



az, x r x p(V + azx r)dV^ 



+ 



Is, 



El 

8tt 



(r x n) z dS A 



[r x (az x r)] zP V • fidSU, (14) 

fs A 

where V = u — az x r is the velocity vector in the frame 
rotating with angular velocity az, z is the unit vector that 
points in the ^-direction, Sa is the Alfven surface that de- 
limits the volume Va , and n is the normal unit vector to the 
Alfven surfac^E The first term on the right of Eq. ([T4|) does 
not contribute to the ^-component torque and is therefore 
null. The second term disappears in the case of a spheri- 
cal Alfven surface. It is also null in the axi- symmetric cases 
considered here (i.e., the dipolar cases), but it is non-null in 
the cases where a surface magnetic map is considered and it 
becomes relatively more important for the cases with larger 
adopted /3o. The third term is the dominant term in our 
simulations. 

We can estimate the time-scale for rotational braking 
as r = J I J, where J is the angular momentum of the star. 
If we assume a spherical star with a uniform density, then 
J — 2/5M*Rl£lo and the time-scale is 



9 x 10 3 



J 



Mo 



Myr, 



(15) 



Po) Uo, 

where Po = 2tt/Qo is the rotational period of the star. For 
V374 Peg, this results in 

6.45 x 10 35 _ 

r ~ = Myr. (16) 



J 

Because J depends on the mass flux crossing a given 
surface, i.e., on the mass- loss rate of the wind M, from 
Eq. ([T3]) . we have a rough scaling relation between J and 
M for cases IDip/lMap, 2Dip/2Map, and 3Dip/3Map 



J oc M oc n. 



1/2 



(17) 



1 The computation of the local unit vector n normal to the Alfven 
surface was done w ith the routine developed by F. Alouani Bibi 
dQpher et alJl2009h . 
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which implies in a time-scale [Eq. (|16[)] for rotational braking 
that scales as 

r(xn~ 1/2 . (18) 

For cases IDip, 2Dip, and 3Dip, this implies that r c± 
28n~ 1/2 Myr, while for cases IMap, 2Map, and 3Map, 
r ~ 18n 12 7 Myr. In either case, r is we ll below the es- 
timat ed solar spin-down time tq c± 7 Gyr (| Weber &; Davisl 

Table [T] presents the mass and angular momentum loss 
rates, and the time-scale for rotational braking calculated 
for all simulations, where we verify the approximate scal- 
ing given by Eqs. (fT3]h JT7}, and JT5J). Comparing to the 
observationally derived rotational braking time-scales of a 
couple of hundreds of M yr for dM stars is open clusters 
(jScholz fe E isloffel 2007), we tend to rule out cases with 
larger coronal base densities (i.e., no > 10 11 cm -3 ). Ac- 
cording to this comparison, the most plausible wind density 
is the one used for models IDip/lMap. Such a density is 
also able to reproduce typical emission measures of dM stars 
(EM & 10 51 cm -3 ) and comparatively (with the remaining 
cases) smaller mass-loss rates and higher wind velocities. 

Because of its high mass-loss rate, V374 Peg is expected 
to lose about 0.01 M (i.e., ~ 3.5% of its mass) during a 
time-scale r. If the star persists on losing this amount of 
mass during a longer period of time, this could have signif- 
icant effects on its evolution. However, when the star ages 
and both its rotation as well as the stellar surface magnetic 
field intensity are reduced, its mass- loss rate is expected to 
become less intense and therefore present a less significant 
effect on the evolution of the star. 

There is also great interest in low-mass stars at 
young ages, when accretion discs are still observed. Accre- 
tion rates observed for pre - main sequence low mass stars 
(jHerczeg fe Hill enbrand 120081 ) are comparable to the wind 
rates that we derived here. Although it is not straightfor- 
ward to extend our models to these very early phases, we 
note that the high wind rates we predict could have impor- 
tant consequences for the evolution of low-mass stars and 
their discs. 

With the inclusion of an observed distribution of sur- 
face magnetic field, we head towards a more realistic mod- 
elling of magnetised coronal winds. Never the less, our model 
presents limitations, such as the neglect of a detailed en- 
ergy balance. Instead, we consider a polytropic relation be- 
tween pressure and density (or temperature) parametrised 
through 7 in the derivation of the energy equation Eq. Q. 
Once the magnetic field distribution is set, the thermal pres- 
sure adjusts itself in order to provide a distribution of heat- 
ing/cooli ng that is able to sup port the MHD solution ob- 
tained (|Low fe Tsinganosl[T986h . If the wind of V374 Peg is 
able to cool down, e.g., by radiative cooling not considered 
in our models, the terminal velocities of the wind could be 
considerably smaller. Depending on where in the wind en- 
ergy deposition (or removal) occurs, the wind velocity may 
change, without affecting the mass-loss rates. For instance, 
if a substantial cooling occurs above the Alfven surface, the 
velocity profile of the wind from that point outwards will 
be affected. As the information of what is happening above 
the Alfven point cannot be transmitted to the sub-Alfvenic 
region, the wind density and velocity profiles in the proxim- 




Figure 4. Isosurface of \B r \/B = 0.97 calculated for the MHD 
model (case 3Map). Color represents the radial velocity of the 
wind. 

ity of the star will not be changed, and consequently neither 
the stellar mass- loss/ angular momentum- loss rates. 

4.2 Departure from Potential Field 

Here, we compare the results of our MHD modelling of the 
wind of V374 Peg with the output of the PFSS method, 
which was used as the initial configuration for the magnetic 
field lines in Set 2 of our simulations. Such a configuration 
is shown in Fig. [T}). The spherical source surface is located 
at r = rss = 5 R* and represented by the dashed circum- 
ference. The field is considered to be everywhere potential, 
and its surface distribution derived from observed surface 
magnetic maps. Above r = rss, the PFSS method assumes 
that the magnetic field is purely radial, decaying as r -2 . Ide- 
ally, we would like to compare the size of the source surface 
derived from the MHD modelling. However, this compar- 
ison is not straightforward, as there is no surface beyond 
which the magnetic field lines are purely radial in the MHD 
model. To overcome this difficulty, Figure |4] shows an isosur- 
face of | B r \/B = 0.97, which represents the surface where 
97 percent of the magnetic fi eld is contained in the radial 
component (|Rilev et al.ll2006l ). The first immediate conclu- 
sion from the comparison between this isosurface and the 
source surface of the PFSS method is the lack of sphericity 
of the former. Moreover, the size of the surface is consid- 
erably smaller than the radius rss = 5 i?* adopted in the 
PFSS method. In our MHD model, the average size of the 
isosurface is (rss) = 2.6 i?*, extending toward the poles out 
to 1.9 i?* and toward the equator out to 2.9 i?*. We note 
that this size is very similar to the adopted size of 2.5 solar 
radii from PFSS methods of the solar coronal magnetic field. 

The stored magnetic energy contained in the potential 
field is in the lowest state, i.e., it is the the minimum value of 
energy that the magnetic field lines can store. In the MHD 
wind case, excess energy is contained in the magnetic field 
lines due to stresses imposed by the wind. To quantify the 
departure of the MHD solution from the potential field so- 
lution, we evaluate the stored magnetic energy in each case. 
Defining the mean magnetic energy as 



10 A. A. Vidotto et al. 




Figure 5. Current density. 



(B 2 ) 



J v B 2 dV 



(19) 



where V is a given volume, we calculated the ratio / between 
the energy contained in the MHD solution and the one in 
the PFSS solution 



/ 



{B 2 h 



(B 2 h 



(20) 



For a volume contained between r = 1 and 5 i?*, / = 1.5, 
implying that the MHD solution has stored about 50 percent 
more magnetic energy than the potential field solution. The 
same evaluation from r = 1 to 2 R* results in / = 1.3 and 
from r = 2 to 5 i?* results in / = 2.8. We conclude that 
closer to the star, the MHD solution deviates little from 
the potential field solution. However, the departure from a 
potential field becomes more important farther out from the 
star. 

The existence of currents in the MHD solution also il- 
lustrates the departure of the magnetic field configuration 
from a potential one. In contrast to the PFSS method, which 
assumes null currents everywhere (j oc V x B = 0), in the 
MHD wind modelling, current densities arise naturally. Fig- 
ure [5] shows the profile of current densities, which is similar 
in all the simulations of the wind of V374 Peg. This simi- 
larity occurs because, in the ideal MHD, j only depends on 
the intensity and configuration of the magnetic field, whose 
characteristics are similar in our simulations. The current 
density is maximum near the star and it decreases with dis- 
tance. A more intense current sheet is seen separating the 
closed-field line region and the open-field line region. 



4.3 Effects on Habitability 

The concept of habitable zones relies on the location where 
liquid water may be found, which is essentially determined 
by the luminosity of the star. Because dM stars are low- 
luminosity objects, the habitable zone is believed to lie near 
the star. However, a strong wind is able to erode the at- 
mosphere of a planet, removing an important shield to the 



creation and development of life. The ram pressure p w of the 
wind that impacts on a planet orbiting V374 Peg is 



pu r 



mu r r 



(21) 



where m — pu r r . In the asymptotic limit of the wind ve- 
locity, for large distances from the star, 

• " 2 (22) 



4 7rr 2 



Using our scaling relations for V374 Peg (cases IDip/lMap, 



2Dip/2Map, 3Dip/3Map) M ~ 4 x 10 



-10 1/2 



M 



yr 



^oo ~ (1500 - 2300)r 



-1/2 



Pw, oo 



(3 to 4.6) x 10 2 



km s \ we have 
dyn cm -2 , 



_1 and 



(23) 



where the range of values exists due to different latitudi- 
nal wind velocities. We note that p w ,oo does not depend 
on the density assumed for our wind models of V374 Peg. 
Therefore, even an increase in the wind mass- loss rate (e.g., 
through an increase in the wind base density) would not 
have any further effect on an orbiting planet, because there 
is an associated decrease in the wind velocity. 
For the solar wind 

^18 



Pw ,oo ,0 ' 



QUoo,Q 

4 7rr 2 



4 x 10 1 



dyn cm 



(24) 



where we have used Mq = 2 x 1O _14 M0 yr -1 and Uoo,q — 
400 km s _1 . Therefore, the wind of V374 Peg has a ram 
pressure that is about 5 orders of magnitude larger than for 
the solar wind. The immediate question we ask ourselves is 
how strong the planetary magnetic field should be in other 
to shield the planet's atmosphere from the erosive effects of 
the wind. Pressure balance between the wind total pressure 
and the planet total pressure requires that, at a distance tm 
from the planet, 



p w B w 



47T 



+ Pp, 



(25) 



where p w , B w and p are the local wind ram pressure, mag- 
netic field and thermal pressure, and B v , rM and p p are the 
local planetary magnetic field intensity and thermal pres- 



sure. Far from the central star, p u 



Pu 



while B w and 



p become negligible. Assuming that the planetary magnetic 
field is dipolar, B P ^ M — B p (R P /rM) 3 , where B p is the equa- 
torial planetary magnetic field intensity at its surface and R p 
is the planetary radius. A minimum magnetic shielding of 
the planet requires that tm > R P . Therefore, for a planet 
orbiting V374 Peg at r = 1 AU, the minimum planetary 
magnetic field required is 
1/2 



B 



(27T p u 



0.1 G, 



(26) 



where we assumed that p p is negligible. A significant p p 
would only increase the limit obtained before. Therefore, 
for a planet orbiting far from the star (~ 1 AU), a mag- 
netic field of only 0.1 G would be sufficient to prevent 
atmospheric erosion caused by the wind. In the asymp- 
totic limit of the wind, Equations ([22]) and ([26]) show that 
£p,min — 0.1[(1 AU)/r] G. 

However, it should be noted that for planets orbiting 
closer to the star, p w > p w ,oo and both the magnetic and 
thermal wind pressures in Equation ([25]) may not be negli- 
gible. Therefore, the minimum planetary magnetic field for 
shielding is 
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£ P ,min = (2irp w + B 2 W + 4n(p - p P )) 1/2 . (27) 

The habitable zone for a 0.5M© star is believed to lie be- 
tween 0.2 - 0.5 AU, while for a 0.1M© star , it is thought 
to lie much closer, between 0.02 - 0.05 AU (Kas ting et al.l 
1993). Therefore, using the results from our wind models 
and neglecting p p , for a planet orbiting on the equatorial 
plane of V374 Peg at r = 0.03 AU, 

5p,min — 4 G. (28) 

This implies that, in order to prevent atmospheric erosion, a 
planet orbiting inside the habitable zone of V374 Peg would 
require a minimum magnetic field intensity of a few G. Con- 
sidering that this is the magnetic field at the equator of the 
planet, its minimum field intensity at the pole is ~ 8 G, 
which is abou t half the intensity o f Jupiter's magnetic field 
(~ 16 G, IConnernev et al . 1998). In fact, a planet with 
the same magnetic field intensity as Jupiter would have, 
at 0.03 AU from V374 Peg, a magnetospheric radius of 
tm — 1.5 R P , where we assumed that p p ~ 0. 

Ultimately, when the star ages, the stellar rotation 
brakes, reducing the stellar surface magnetic field intensity, 
and therefore the wind velocity. In this case, the ram pres- 
sure felt by the planet will tend to decrease. 



5 CONCLUSIONS 

We have presented 3D MHD numerical simulations of the 
wind of the M dwarf star V374 Peg. Starting from an initial 
magnetic field configuration and a thermally-driven wind, 
the wind is evolved in time, resulting in a self-consistent in- 
teraction of the wind and the magnetic field. We analysed 
cases where the surface distribution of magnetic field was 
set to be either simply dipolar or set as potential field ex- 
trapolations from observed magnetic field maps. This is the 
first time that an observationally derived surface magnetic 
field map has been implemented in a MHD model of stel- 
lar wind for a low-mass star. Our first goal was to confirm 
that observed maps could be reliably incorporated into the 
MHD code by comparing the results with those from a sim- 
ple dipole. This is a necessary first st ep before the more 
complex fields of other low-mass star s (|Morin et al . 2008a; 
iMorin et al1l2010l : iDonati et alll2008l ) can be modelled. In 
all our simulations, we found MHD solutions for the wind 
of V374 Peg, showing that it is possible to develop coronal 
wind models with a realistic distribution of magnetic field. 
Our solutions showed that, close to the surface of the star, 
the magnetic field deviates little from the potential field so- 
lution. The 'effective source surface' of the MHD models, 
i.e., where the magnetic field is mainly cointained in the ra- 
dial component, has an average radius of 2.6 itU, about half 
the size of the radius adopted in the potential field solution 
used as initial condition. 

For V374 Peg, the model parameter that is better con- 
strained is the surface magnetic field, from which magnetic 
maps have been derived. Other parameters needed by our 
models, such as coronal densities and temperatures, were 
adopted as representative of dM stars, and varied as to anal- 
yse their implications on the wind of V374 Peg. Adopting 
such parameters (Table[T]), we noted that V374 Peg presents 
a highly magnetised corona, with plasma- j3 parameters of 
the order of 10 -5 - 10 -3 at the base of the corona. This is 



considerably smaller than f3 adopted for the base of the solar 
corona (/3© ~ 1), which implies that the wind of V374 Peg 
may deviate greatly from a low-velocity, low-mass-loss rate 
solar- type wind. In the fast magnetic rotato r limit, the more 
simplistic stellar wind model developed by IWeber fe Davisl 
dl967h predicts terminal velocities of ~ 3320 km s 1 for a 
wind mass-loss rate of about 10 -11 M© yr _1 ( Appendix [A]) . 

For our models with T = 2 x 10 6 K and B = 1660 G, 
we showed that terminal velocities are around (1500 - 
2300)n^ 2 1 ^ 2 km s _1 where ni2 = no/(10 12 cm -3 ), while the 
mass- loss rates are about M ~ 4 x 10" 10 n^ 2 M yr" 1 . We 
also estimated the characteristic time r that the magnetised 
wind of V374 Peg would take to remove angular momentum 
of the star, and therefore, brake stellar rotation. We found 
that the angular momentum-loss rate scales with the mass- 
loss rate, and that r ~ 28n^ 2 1//2 Myr. Compared to observa- 
tionally derived values from period distributions of stars in 
open clusters, this value suggests that V374 Peg may have 
low coronal base densities (< 10 11 cm -3 ). Once wind den- 
sities are better observationally constrained, we will be able 
to further restrict our model parameters and, consequently, 
the derived wind characteristics. The scaling relations we 
found should be applicable to other highly magnetized dM 
stars, although their constants should differ from the ones 
obtained for V374 Peg. 

We investigated the possible effects the wind may have 
on the habitability of a planet orbiting V374 Peg. We found 
that the wind ram pressure impacting on a planet is about 
5 orders of magnitude larger than the ram pressure of the 
solar wind for a planet orbiting at same distance. This high 
ram pressure of the wind of V374 Peg could result in ero- 
sion of the planet's atmosphere in the case the planet does 
not possess a sufficiently high magnetic field able to shield 
it from the stellar wind. For a planet orbiting at ^ 1 AU, 
we estimated the minimum necessary magnetic field to be 
~ 0.1 G. However, for planets much closer to the star, pos- 
sibly in the habitable zone, the wind ram pressure is much 
larger, which could present a problem for the creation and 
development of life. Never the less, we estimated that a plan- 
etary magnetic field of about half Jupiter's intensity would 
be able to protect the planet from the erosive effects of the 
powerful stellar wind of V374 Peg. 
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APPENDIX A: COMPARISON TO SIMPLER 
WIND MODELS 

It is interesting to have in mind the wind characteristics 
expected from more simplistic wind mo dels. Non- rota ting, 
non-magnetised, thermally-driven winds (|Parkerlll958h pre- 
dict ter minal velocities that are of t he order of the Parker 
velocity (|Belcher fe MacGregorlll976h 

/ ?r 2 \ 1/2 

parker / ^ u s,0 2 , 2 \ /A1\ 

u ™ = Vt-^l Vesc, ° Ur, ° J ' 
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where c s ,o = [yks To/ ({imp)] 1 / 2 is the sound speed, v esc ,o = 
(2GM*/i?*) 1/2 is the surface escape velocity, and u r ,o is 
the radial wind velocity at i?*. Neglecting u r ,o in face of 
the other velocities in Eq. (|Alj) . for V374 Peg, u^ rker ~ 
290 km s" 1 for T = 2 x 10 6 K, and u^ rker ~ 1290 km s" 1 
for T = 10 7 K. 

Because V374 Peg is a highly magnetised star, spinning 
with a rotational period of about ~ 11 h, V374 Peg is a fast 
magnetic rotator, where magneto- cent rifugally force pro- 
vides an efficient mechanism to accelerate t he wind. On the 
basis o f the stellar wind model developed by I Weber &; Davis] 
(1967) for a magnetised (with a radial magnetic field), ro- 
tating star, the wind terminal velocity in th e fast magnet ic 
rotator limit is given by the Michel velocity (|Michel|[l969h 



michel 



Qo(^r,oi?*) 2 1 



(A2) 

M J 

In such a limit, u™ chel > c s ,o and u™ chel > u£ rker . For 
V374 Peg, 

^michel _ 3320M-V3 km s -l } ( A3 ) 

where Mn = M/(10~ n M yr _1 ). For cases IDip/lMap, 
where Mn « 4, u™ chel ~ 2100 km s" 1 , while for cases 
3Dip/3Map, where Mn « 40, <f hel - 970 km s _1 . It 
is interesting to compare these values with the ones ob- 
tained in our simulations for cases IDip/lMap, 2Dip/2Map, 
and 3Dip/3Map (i.e., Uoo « [1500 - 2300]n~ 1/2 km s" 1 ). 
Compared to the corresponding wind terminal velocity Uoo , 
^michel j g a ]3 OU ^ 7 to 11 times smaller than Uoo for cases 
IDip/lMap, and about 1.6 to 2.4 times smaller for cases 
3Dip/3Map. 
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